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io The rapidity dependence of two-particle momentum correlations can be used to probe the viscosity 

ii of the liquid produced in heavy nuclei collisions at RHIC. We reexamine this probe in light of 

12 the recent experimental analyses of the azimuthal-angle dependence of number correlations, which 

o demonstrate the importance of initial state fluctuations propagated by hydrodynamic flow in these 

14 correlations. The NEXSPHERIO model combines fluctuating initial conditions with viscosity-free 

O 1 15 hydrodynamic evolution and, indeed, has been shown to describe azimuthal correlations. We use 

D i6 this model to compute the number density correlation 7?2 and the momentum current correlation 

^J 17 function C, at low transverse momentum in Au+Au collisions at ^/snn = 200 GeV. The correlation 

_^> i8 function C is sensitive to details of the collision dynamics. Its longitudinal width is expected to 

I 19 broaden under the influence of viscous effects and narrow in the presence of sizable radial flow. 

20 While the NEXSPHERIO model qualitatively describes the emergence of a near-side ridge-like 

i— i 21 structure for both the R2 and C observables, we find that it predicts a longitudinal narrowing of the 

wG 22 near side peak of these correlation functions for increasing number of participants in contrast with 

recent observations by the STAR Collaboration of a significant broadening in most central collisions 

relative to peripheral collisions. 
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(N 25 I. INTRODUCTION 

> 

qq 26 Measurements of elliptic flow indicate that strong collective flow is achieved in nuclear collisions at RHIC [IH5]; 
l/"~) 27 see also the reviews [6j [7]. Hydrodynamic models describe many features of the experimental data, indicating that 
C*") 28 the produced medium behaves as a nearly ideal fluid, i.e. a fluid endowed with a very small shear viscosity per unit 
j^| 29 of entropy, rj/s. Considerable effort is devoted to using elliptic flow data to measure 77/s; see, e.g., [El IS], and Refs. 
f^*) 30 therein. However, the interpretation of such measurements has many sources of theoretical uncertainty. The primary 
r—i 31 uncertainties are: initial conditions, event-by-event fluctuations, and freeze-out |10j . These effects are somewhat 
"^ 32 outside the scope of hydrodynamics and further theoretical work is required to elaborate a systematic scheme to 
^* 33 compute them with quantifiable uncertainties. It is therefore important to explore independent methods to determine 
77/s. 

One such alternative method was proposed by Gavin et al. based on two-particle transverse momentum correlations 
[llj . This novel technique relies on the fact that viscous forces impart momentum transfers between different fluid cells. 
These momentum transfers in turn modify the relative rapidity of correlated particles produced by these cells thereby 
leading to a longitudinal broadening of transverse momentum two-particle correlation functions. The longitudinal 
profile of this correlation function, hereafter called C, is expected to broaden with increasing collision system life- 
time, and for increasing number of collision participants. The measurement strategy proposed by authors of [llj is 
to measure the px correlation function C as a function of collision centrality and seek evidence for its longitudinal 
broadening in most central collisions. Given its explicit dependence on particle momenta, the correlation function C 
is sensitive to momentum current correlations and as such constitutes a more suitable observable to measure viscosity 
than simple two particle (number) correlation functions. 

The STAR experiment recently carried an analysis of a differential version of the pr correlation function C. STAR 
reported that the near-side of the correlation function broadens progressively with increasing collision centrality [12j . 
The broadening is considerable but translates nonetheless in a rather modest rj/s consistent with values deduced from 
comparison of hydrodynamic calculations with vi measurements. While, this agreement is encouraging, one must 
however exercise caution in the interpretation of the broadening of the measured correlation function. Indeed, one 
must assess whether other mechanisms could produce or perhaps mask broadening. Reaction mechanisms and effects 



of particular interest include two-particle kinematical focusing associated with radial flow, jet production, and more 
specifically jet quenching, flux tubes, fluctuations, anisotropic flow, resonances, and momentum conservation. One 
must also acknowledge that the structure of two-particle correlations observed is not fully elucidated. While differential 
correlation functions have been previously measured by the STAR |13|. PHENIX [T3] and PHOBOS [T5] experiments, 
a significant interest lies in studying diverse correlation functions which provide better understanding and sensitivity 
to various reaction mechanisms. Dihadron correlations have previously been studied containing a high-pT "trigger" 
particle and a low-p^ "associated" particle. Measurements of near-side enhanced yield - long range correlations - also 
known as the "Ridge" were reported by the STAR collaboration [THl Q~7] . Some theoretical models attribute the ridge 
to jet-medium interactions: particles from jet fragmentations in QCD vacuum result in a peak at At? ~ and those 
affected by the medium are diffused broadly in pseudorapidity, forming a ridge-like structure [18-25 . Other plausible 
explanations of the ridge formation are given by transverse radial flow models where kinematic focussing of clusters, 
strings, or color flux tubes is induced by increasing radial flow with centrality |26j . Transverse radial expansion 
produces strong position-momentum correlations that leads to characteristic rapidity, px and azimuthal correlations 
among the produced particles. Another interesting and promising interpretation of the near-side ridge observed in both 
triggered and untriggered correlation functions is based on the notion of initial density fluctuations. Initial density 
fluctuations in the coordinate space shall lead to fluctuations in the number of locally produced particles. These 
fluctuations translate, under the influence of hydrodynamical radial expansion, into long range pseudorapidity two- 
particle correlations. Such long range correlations are well demonstrated in the NEXSPHERIO simulations reported 
by Takahashi et al. 27J. Other studies based on hydrodynamical and hybrid models have essentially confirmed this 
important observation J25H52] , Also note that much work is being done to quantatively understand the connection 
between initial collision system spatial anisotropies and final state flow coefficients, v n , of orders n = I through 6 
[55H57] . and to determine the fluid viscosity (see for instance [3J3I3S])- NEXSPHERIO describes many features of 
the measured correlations, particularly their azimuthal dependence. It is thus of interest to extend the work of these 
authors and consider predictions of NEXSPHERIO for transverse momentum two-particle correlations. 

The NEXSPHERIO model, the correlation observable of interest, and the analysis technique used in this work are 
described in the next section. Results are presented in sec. |III| and conclusions summarized in sec. |IV| 



27 II. NEXSPHERIO MODEL AND ANALYSIS TECHNIQUE 

28 NEXSPHERIO gDl S5 combines the event generator NEXUS g2] and the hydrodynamical code SPHERIO 03] to 

29 simulate particle production in A + A collisions. NEXUS provides for initial parton-parton interactions and energy 

30 deposition with lumpy initial conditions (ICs) specified in terms of number density, energy density and baryon density. 

31 SPHERIO is then used to simulate the hydrodynamical evolution in each NEXUS event, i.e thereby including the 

32 effect of IC fluctuations in the production of collective flow patterns and correlated particle production. Indeed, the 

33 recent studies with NEXSPHERIO show that hydrodynamic expansion of fluctuating ICs lead to the formation of 

34 near-side ridge-like structure in At? vs. A<p correlation functions |27j while, by contrast, smooth initial conditions 

35 do not feature any near-side ridge-like structures. Furthermore, an interesting interpretation of the origin of the 

36 away-side double peaked structure of the two particle correlation function was also discussed in the context of the 

37 NEXSPHERIO lumpy initial condition in reference [H] . 

38 The transverse momentum two-particle correlation C depends at a basic level on particle number fluctuations and 

39 should therefore exhibit the formation of a ridge, in NEXSPHERIO, for central collisions with lumpy initial conditions. 
« C also depends on actual momentum-momentum correlations. Its magnitude and dependence on Ai] vs. Acj> should 
4i thus also be influenced by radial hydrodynamical flow brought about by lumpy initial density conditions. However, 

42 given SPHERIO is built as an ideal hydrodynamical transport code, there should be no viscous effects: no broadening 

43 of the correlation function should arise with increasing collision centrality. Any change in the correlation function 

44 profile and width should be due solely to viscous free expansion dynamics. 

45 Our analysis is based on a total of 500,000 simulated Au + Au collisions events at y'sjvjv = 200 GeV produced with 

46 NEXSPHERIO [27]. Events were grouped in four collision centrality classes corresponding to 0-10% (most central), 20- 

47 30%, 40-60% and 60-80% (most peripheral) fractions of the total reaction cross-section. Modeling nuclear collisions in 

48 ideal hydrodynamics is limited to low pt (< 2.0 GeV/c) because it is generally assumed that the equilibrium description 

49 fails at high px, where particle production may be dominated by hard processes. We therefore focus our analysis 
so on bulk particles, i.e. particles produced in the range px < 2.0 GcV/c. This momentum selection is used for both 

51 particles in constructing the two particle correlation function i?2 , and the momentum dependent correlation function, 

52 C, defined below, thereby yielding "untriggered" distributions. We simulate the STAR experiment acceptance and 

53 further limit our simulations to particles produced in the pseudorapidity range |?7| < 1.0 and with 0.2 < px < 2.0 

54 GeV/c. 

55 We carried an analysis of observables i?2 and C studied as a function of pseudorapidity (Arf) and azimuth(A</>) 



differences of particle pairs. These correlation functions are defined as follows: 
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where p\ — dNi/dpTdrjd(j) and p 2 = dN 2 j dp t ^\dr\\dfadpt^ 2 dr\ 2 d4) 2 correspond to single and two-particle densities, 
respectively; {pt) k = C^2pt,i)u I ( n )k * s the average particle momentum. The label k=l,2 stands for particles being 
detected at (pt./c^fcj^fc), and the brackets represent event ensemble averages. The integrals are taken over 2tt az- 
imuth and in the same pseudorapidity range for both particles. The delta functions are used to restrict correlation 
yield contributions to relative azimuthal angle and pseudorapidity differences only. The observable R 2 measures the 
correlation between the number of particles emitted at relative pseudorapidity Arj and azimuthal angle difference A<j). 
It is a robust observable and is expected to scale inversely to the number of correlated particle sources in the absence 
of collective effects and rescattering of secondaries. The observable C amounts to a differential version of the integral 
correlation used by authors of [11] . It is sensitive not only to the degree of correlation between produced particles but 
to the hardness of particle spectra, and correlations between the momentum of particles. 

III. RESULTS 

We present, in Fig. [II the correlation function R 2 for (a) 60-80% (most peripheral), (b) 40-60%, (c) 20-30% and 
(d) 0-10% (most central) collision centralities. Our calculation of R 2 reproduces the features already reported in [STJ 
for a 'per trigger' correlation function. R 2 exhibits a somewhat narrow near-side peak (i.e. near Arj f=a A</> ss 0) and 
an extended ridge-like structure centered at A<p ss n in all collision centralities. The correlation function exhibits a 
strong cos{2Afa modulation associated with collective flow in mid-central collisions. This component is minimal in 
most central collisions but R 2 exhibits the build up of a narrow ridge at A(f> ~ which extends over the full Arj range 
of the correlation function. 

Figure [2] presents the correlation function C for the same four collision centrality bins used for R 2 in Fig. [I] C 
exhibits correlation shape dependencies on A 77 and A<fi qualitatively similar to those observed for R 2 . It also has a 
narrow near-side peak at Arj f=a A(f> rj and an extended Arj ridge at A<f> ~ 7r in all collision centralities. Mid-central 
bins feature strong elliptic flow modulations and the most central bin exhibit a near-side peak atop an extended ridge. 
Also note that both R 2 and C exhibit an overall decrease of the correlation function strength consistent with the 
increase in the average number of participant nucleons and the decrease of the elliptic flow modulation from 60-80% 
to 0-5% collision centrality ranges. 



(a) 60-80% 




FIG. 1: (Color online) Two-particle correlation function, R2, plotted as a function of the particle pair pseudorapidity difference, 
A??, and relative azimuthal angle, A(j> [rad], for (a) 60-80%, (b) 40-60%, (c) 20-30% and (d) 0-10% collision centralities in Au 
+ Au interactions at ^smn = 200 GeV calculated with the NEXSPHERIO event generator. 



29 Authors of Ref. [27 showed the initial density (event-by-event) fluctuations present in NEXSPHERIO calculations 

30 produce a near-side ridge in 'triggered correlation functions'. We here find that these initial density fluctuations also 

31 produce a near-side ridge in both 'untriggered' and transverse momentum correlation functions. Given C has an 
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FIG. 2: (Color online) Two-particle pr correlation function, C, plotted as a function of the particle pair pseudorapidity 
difference, Ary, and relative azimuthal angle, A(j> [rad], for (a) 60-80%, (b) 40-60%, (c) 20-30% and (d) 0-10% collision centralities 
in Au + Au interactions at ^/snn = 200 GeV calculated with the NEXSPHERIO event generator. 



explicit dependence on the particle momenta, one expects measurements of this variable should provide additional 
constraints for theoretical models. We further observe that the away-side ridge-like feature evolves significantly from 
a broad structure at Arj ~ in most peripheral collisions to a saddle-like shape with a minimum located at Ar/ ~ 
in the most central collisions. 

The observed build up of a near-side ridge in most central collisions is understood in NEXSPHERIO simulations 
to arise from hot spots extending in the longitudinal (i.e. beam) direction. Fluctuations in number and transverse 
position of these hot spots generate finite azimuthal anisotropy which one expects to be essentially independent of 
pseudorapidity in the acceptance considered in this work. Given that NEXSPHERIO involves ideal hydrodynamical 
transport of the energy deposited during the early moments of the collisions, one could thus anticipate that the 
longitudinal width of the near-side peak should be essentially unaffected by the production of the near-side ridge. We 
verify this assertion by plotting Ar) projections of the R2 and C correlation functions in the range \Acf>\ < 1.0 radians 
in Figures 3 and 4 respectively. We parameterize and fit these distributions using a simple two-component ansatz to 
study their evolution with collision centrality. 
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FIG. 3: Projections of the correlation function R 2 for \A<f>\ < 1.0 radians on the Ar] axis for (a) 60-80%, (b) 40-60%, (c) 20-30% 
and (d) 0-10% centralities in Au + Au collisions at -y/sjvjv = 200 GeV from NEXSPHERIO event generator. 



Fit results are presented in Table I and II for observable R2 and C respectively. We first observe that the near-side 
peak amplitude of both observables decreases monotonically with increasing collision centrality. We also observe the 
emergences of strong azimuthal modulation in mid-central collisions, and the rise of a near-side ridge extending over a 
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FIG. 4: Projections of the correlation function Cfor \A<j>\ < 1.0 radians on the An axis for (a) 60-80%, (b) 40-60%, (c) 20-30% 
and (d) 0-10% centralities in Au + Au collisions at ^/s NN = 200 GeV from NEXSPHERIO event generator. 

long At/ range in the most central collisions. These result from the formation and evolution of the medium, with initial 
fluctuations, and the effects of radial and elliptic flow, embodied in the NEXSPHERIO calculations. We furthermore 
find that the distribution widths, a, exhibit a decreasing trend with increasing collision centrality. Specifically, we 
find that the At/ width of C monotonically evolves from 0.72 ± 0.02 units of rapidity in peripheral collisions, to 0.54 
± 0.03 in most central collisions. R2 exhibits a similar behavior, with an initial width of 0.81± 0.02 in peripheral 
collisions, and reduces to 0.63 ± 0.01 in 20-30% collisions. It however rises slightly in most central collisions. We 
also characterized the width of the distributions by calculating their root-mean-square (RMS) above a constant offset 
evaluated based on the value of the correlation function at At/ ~ ±2. The deduced RMS values are listed in Table I 
and II. They are found to be in agreement, within statistical errors, with the widths, a, obtained from the fits. 

TABLE I: Various fit parameters for R 2 

Centrality Offset Amplitude a RMS X 2 / nd f 

60-80% 0.345 ± 0.001 0.39 0.81 ± 0.02 0.78 ± 0.04 38.32/28 

40-60% 0.156 ± 0.001 0.17 0.69 ± 0.01 0.68 ± 0.03 41.16/28 

20-30% 0.032 ± 0.002 0.038 0.63 ± 0.01 0.63 ± 0.03 24.27/28 

0-10% 0.022 ± 0.002 0.025 0.74 ± 0.03 0.75 ± 0.04 37.79/28 



TABLE II: Various fit parameters for C 



Centrality Offset 

60-80% 0.1096 ± 0.0002 
40-60% 0.0559 ± 0.0001 
20-30% 0.0142 ± 0.0001 
0-10% 0.0084 ± 0.0001 



Amplitude a RMS xV nd f 

0.126 0.72 ± 0.02 0.71 ± 0.04 23.61/28 

0.063 0.67 ± 0.01 0.66 ± 0.03 16.23/28 

0.017 0.56 ± 0.01 0.59 ± 0.03 17.55/28 

0.0097 0.54 ± 0.03 0.58 ± 0.02 33.99/28 



Figure 5] shows a comparison of the RMS obtained with NEXSPHERIO (black squares) with those reported by 
STAR (black circles) for Au + Au collisions at y/s^j = 200 GeV. The RMS predicted by NEXSPHERIO substan- 
tially exceeds the value measured by STAR in most peripheral collisions. We find additionally that NEXSPHERIO 
predicts a progressive narrowing of the correlation function for increasing number of participants, at variance with 
the observation, by the STAR experiment, of an increasing broadening for larger number of participants |12j . The 
RMS reported by STAR for most central collisions is in fact almost twice as large as that predicted by the model. In 
addition, the smaller value of RMS observed in the STAR peripheral data could in part be due to contributions from 
jet particles that will yield a narrower correlation function. 

We thus find that while NEXSPHERIO qualitatively reproduces the azimuthal modulation and near-side ridge 
observed in experimental data it predicts a narrowing rather than the observed longitudinal broadening of the cor- 
relation function C. The longitudinal narrowing in NEXSPHERIO is in large part due to the strong transverse flow. 




part 



FIG. 5: RMS of the correlation function C predicted by NEXSPHERIO (solid squares) as a function of the number of participant 
nucleons, compared with the results reported by STAR [12] (solid circles), for nine centrality classes in Au+Au collisions at 
^s^ = 200 GeV. 



This narrowing effect of transverse flow on rapidity correlations was first discussed in [¥7, 48 . The narrowing found 
here is rather modest, however, amounting to approximately 0.13 units of rapidity. 

To understand why NEXSPHERIO fails to describe the measured broadening, we ask how broadening may arise 
in general. In Ref. [TTj, it was predicted that viscous diffusion can produce a considerable longitudinal broadening 
of the momentum current correlation function C. The idea is that viscous friction acts to reduce the difference in 
transverse flow velocity between neighboring fluid cells. This can be measured by studying the rapidity dependence of 
p t fluctuations, because of the rough correspondence between spatial distance along the beam direction and rapidity. In 
principle, rapidity broadening can also result if the particle production mechanism varies appropriately with centrality. 
We speculate that NEXSPHERIO fails to reproduce the measured broadening because neither effects are included. 

The STAR Collaboration used their measurement of the broadening of the correlation function C to estimate the 
viscosity per unit of entropy based on the following equation [llj : 



T C _S 



'cj 



(3) 



where a c and oq stand for the longitudinal widths of the correlation function in central collisions and at formation 
time, respectively, tq refers to the formation time and t c j is the kinetic freeze-out time at which particles have no 
further interactions |45) . T c stands for a characteristic temperature, here taken to be the critical temperature of 
the medium, T c = 180 MeV. This expression neglects the fact that radial flow causes a narrowing of the correlation 
function C. We estimate the error due to that omission using (3) as follows. A difference in the width due to flow 
increase the value of r\ by a fraction Ary/r/ = Aa 2 /(a 2 — a 2 ), where a CtP are the measured values in central and 
peripheral collisions. NEXSPHERIO yields a reduction of width from 0.72 in peripheral collisions to 0.54 in central 
collisions, while the data increases from 0.54 to 0.94, see Fig. 5. STAR reports a range rj/s = 0.06 — 0.21 in [12] , We 
thus estimate an increase of A77/77 « 0.38 for the upper limit of this range. This suggests that radial flow effects have 
in fact a somewhat limited impact on the width of the correlation function C relative to that of viscous effects and 
can thus be neglected, to first order, in the extraction of the fluid viscosity based on STAR data. 



IV. CONCLUSIONS 



In summary, we presented a study of the centrality dependence of the correlation functions R2 and C in Au + Au 
collisions at ^/s^n — 200 GeV based on the NEXSPHERIO model. We find the two observables exhibit qualitatively 
similar shapes in A77 and A(/>, and dependence on collision centrality. Quantitative differences however arise from 
the explicit dependence of C on particle momenta. We find that both observables exhibit a near-side ridge in central 
collisions owing to event-by-event fluctuations in the initial transverse energy deposition profile. We studied near-side 
projections of the near-side (|A0| < 1.0 radians) of these correlation functions and studied their evolution with collision 
centrality. We found the longitudinal width of C exhibits a small decrease with increasing collision centrality owing to 
radial expansion dynamics. The magnitude of this reduction is rather modest. It is unlikely that this reduction would 
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mask correlation broadening associated with viscous effects in real systems. Indeed, we find that the RMS reduction 
observed in this work is rather modest in comparison with the sizable increase of the longitudinal width reported 
by STAR from their measurement of C in Au + Au collisions at -^/sjvtv = 200 GeV. From this study, it is clear 
that the observed increase in the width of the longitudinal correlation function C from peripheral to central collisions 
by the STAR experiment is not due just to variations of flow effects. In fact, our result shows that with an ideal 
hydrodynamic scenario, the width of the correlation function would decrease, thus suggesting that the broadening 
caused by other possible effects such as viscosity could be even higher than the observed values. Taken at face value, 
the STAR data suggest viscous effects are indeed finite and perhaps well above the quantum bound in central Au + 
Au collisions. 

While NEXSPHERIO does not describe the measured increase in A77 with increasing centrality, the model does 
describe the azimuthal dependence of correlations qualitatively well. We speculate that an extension of NEXSPE- 
RIO that includes viscous hydrodynamics, or other similar viscous hydro-models, might describe the A77 broadening 
reported by STAR. We surmise that this is a generic feature of models that derive long range rapidity correlations 
from the early time behavior followed by hydrodynamic flow [23 I49H52] . We have not dwelled on the detailed ways in 
which measurements of beam-energy, p^, and projectile-mass dependence might ultimately distinguish these models. 

This work was supported in part by U.S. NSF grant PHY-0855369 and DOE grant DE-FG02-92ER40713. We also 
acknowledge support from FAPESP and FAPERJ CNPq, CAPES and PRONEX of Brazil. 
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